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In nonlinear disordered Hamiltonian lattices, where there are no propagating phonons, the spread- 
ing of energy is of subdiffusive nature. Recently, the universality class of the subdiffusive spread- 
ing according to the nonlinear diffusion equation (NDE) has been suggested and checked for one- 
dimensional lattices. Here, we apply this approach to two-dimensional strongly nonlinear lattices 
and find a nice agreement of the scaling predicted from the NDE with the spreading results from 
extensive numerical studies. Moreover, we show that the scaling works also for regular lattices with 
strongly nonlinear coupling, for which the scaling exponent is estimated analytically. This, for the 
first time, shows that the process of chaotic diffusion in such lattices does not require disorder. 



I. INTRODUCTION 

Nonlinearity may have a nontrivial effect on energy 
spreading in extended systems. In regular, homogeneous 
in space media and lattices, nonlinearity can partially 
block or reduce the linear spreading mediated by waves 
(phonons), due to creation of localized structures (soli- 
tons, breathers). Contrary to this, in systems where there 
are no propagating linear waves, nonlinearity is the only 
mechanism responsible for the spreading of initially lo- 
calized wave packets. This situation will be addressed 
in this work. One of the possibilities to suppress the 
linear spreading is to introduce disorder in the system. 
Then, in one and two dimensions, all linear eigenmodes 
become localized and their spectrum becomes purely dis- 
crete^his effect is known as Anderson localization (see [ij 
and 0-01 for recent experimental observations). In this 
case nonlinearity leads to a weak subdiffusive spreading 
of the wave packets, as has been demonstrated in one 0- 
[Tsj and two @, [l^ dimensions. It should be noted that 
these results are mainly based on numerical experiments, 
while purely theoretical attempts (cf. [Tsl - flTj ) have not 
been fully successful in explaining numerical observa- 
tions. Also, by studying chaos properties [isj and the 
possible existence of KAM tori flSl, some reasoning for 
a slowing down of spreading has been proposed. 

Another situation where linear waves are absent are so- 
called strongly nonlinear lattices (20l - [2^ . These lattices 
consist of linear or non-linear oscillators which are cou- 
pled to nearest neighbors by nonlinear forces. The propa- 
gating modes here can be nonlinear waves only, typically 
these waves are compactons [l^. An example of such 
a lattice is the well-known toy "Newton's cradle" p5| . 
Again, disorder in such a lattice blocks the compactons 
(if they exist) and the observed spreading is not a trav- 
eling wave phenomenon, but a slow subdiffusive process. 

All the numerical evidence in the cited literature in- 
dicates for a subdiffusive spreading of wave packets due 
to weak chaos. A natural question is whether this phe- 
nomenon can be described phenomenologically as a cer- 



tain "universality class" , similar to successful statistical 
approaches to problems like percolation and roughen- 
ing of interfaces. In recent papers [H, [l^ such a phe- 
nomenological description based on the properties of the 
Nonlinear Diffusion Equation (NDE) has been proposed 
and tested for one-dimensional lattices. Here, we extend 
this theory to two dimensions, and compare the scaling 
following from the NDE with numerical results from an 
extensive study on 2D strongly nonlinear lattices. We 
emphasize that within this framework we also can inves- 
tigate spreading in regular lattices, where all oscillators 
are in resonance. This is a new, theoretically important 
case, helping to promote the understanding of nonlinear 
spreading, as here also a theoretical prediction on the 
spreading exponent is possible. 

We start in section [ll] by describing the two- 
dimensional Hamiltonian lattices. Then, in section IIIIl 
we introduce the two-dimensional NDE and deduce the 
scaling and spreading properties from its self-similar so- 
lution. This is followed by section IIVI with our main 
results, the comparison of these predictions with the ex- 
tensive numerical simulations including theoretical pre- 
dictions of the spreading exponents when possible. 



II. 2D STRONGLY NONLINEAR LATTICES 

The model considered here is a straightforward gen- 
eralization of the system considered in [l^, HI] to two 
dimensions; it consists of oscillators with power-law lo- 
cal and interaction potentials described by the following 
Hamiltonian; 

H — ^ Ei^k 

i,k 



E.,.k = ^ + — —h.k'^+ 

2 K 

+ ^(Ift+i./c - it.k\^ + \<ii-i,k - qi.k\^ 
+ \qi,k+i - qi,k\^ + \qi,k-i - qi,k\^) 
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1. Instantaneous logarithmic local energy density 
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logj^gWi^k for fi: = 2, A = 4 at times 10^, 10'' and 10" (left 
to right panels). The upper row shows results for a regular 
lattice {oJi^k = 1) with energy E — 1 {W = 13—1 from vari- 
able transformations). The lower row is from simulations of 
a disordered lattice {uJi^k £ [0, 1]) with energy E = 10. The 
total size of the squares is 160x160 lattice sites. 



Here qi^k and pi^k are the positions and momenta of the 
oscillator at lattice site i, fc of a quadratic lattice with 
nearest neighbor coupling. Parameters k and A repre- 
sent the power of the local and the coupling potentials; 
uji^k is a parameter of the local potential and we study 
this model with a random local potential (randomly iid. 
'-^i,k S [0, 1]) as well as in the regular case where uji^k — 1- 
Parameters W and (3 describe the local and the coupling 
strength in this model. 

We have a freedom of rescaling the Hamiltonian and 
the time to get rid of some parameters (this is done sim- 
ilarly to the one-dimensional case [1^ HI]). In the case 
of different nonlinear powers k 7^ A, parameters W and 
/? can be set to = /3 = 1 by rescaling q, p and t, and 
the only remaining parameter is the total energy E in 
the system. 

For a homogeneous nonlinearity with equal powers 
AC = A, by a proper rescaling, the energy and the local 
nonlinear strength can be set to E = W = 1 and (3 is the 
only remaining parameter in this case, describing relation 
of the coupling and the local potential. This possibility 
to scale the total energy to unity induces the following 
scaling relation between the energy and the time [l^, [IJ : 



(2) 



We will use this result later for comparison with the scal- 
ing predictions of NDE. 

Starting from initially localized excitations in such a 
two-dimensional lattice, we find sharply localized spread- 
ing states. This observation is illustrated in Fig. [T]where 
we show the time evolution of a random initial excita- 
tion of 5x5 sites for a lattice with k = 2, A = 4 and 
random local potential {uJi^k G [Ojl]) as well as the reg- 
ular case {uji^k = !)• Our fundamental observable is the 



distribution of energy density Wi^k = Ei^k/E where Ei^k 
is the local energy at site i,k defined in ([T]). One ob- 
serves in Fig. [T] that the extent of the excitation clearly 
increases with time and has very sharp tails at the bound- 
ary (note the logarithmic scale of the grey coding). The 
states spread roughly circular, although with disorder 
(lower row in Fig. [T]) the boundary is more rough than in 
the regular case. We note that, opposed to the ID case 
[20I [231 , no compactons or quasi-compactons that propa- 
gate through the lattice are observed in the 2D case, even 
without disorder (upper row in Fig.[T]). This allows us to 
apply the same approach to both regular and disordered 
lattices. 

To quantify the spreading of energy density shown in 
Fig. [T] we introduce the 2D second moment calculated as: 



An' 



E 

i.k 



{{m^ - iY + {my - kf) Wi^- 



(3) 



with rrix = J2i,k * ^-^"^ "^y = J2i.k ^ ^^,k being the 
center of the distribution. lS.n^ is a measure for the ex- 
citation area in terms of the number of excited sites. In 
order to have a characterization of the uniformity of the 
wave field, we need to calculate other quantities that are 
more or less sensitive to peaks in the distribution. Fol- 
lowing '26] we calculated the Renyi entropies of the distri- 



E»,fc w^,k In w^,k and 5*2 = - In ^, l/< 



bution S - 

(the latter is directly related to the participation num- 
ber), and combined them into the structural entropy de- 
fined as: 

^-str =S-S2. (4) 

If the structural entropy is constant in course of evolu- 
tion, then the relative strength of the peaks in the dis- 
tribution of energy does not change; growing structural 
entropy means that peaks become relatively stronger. 



III. 2D NDE 

The Nonlinear Diffusion Equation (NDE) was found 
to give reasonable predictions on the spreading in one- 
dimensional strongly nonlinear lattices [22.. ,23] . The gen- 
eral picture one has in mind when assuming the NDE to 
describe spreading states is that the motion of already ex- 
cited, active, lattice sites is chaotic. This can be viewed 
as an intrinsic stochasticity of the motion which allows 
for new lattice sites to get excited by stochastic driv- 
ing from the active neighbors. The chaos, and therefore 
also the excitation efficiency, gets weaker in the course of 
spreading, because the energy density decreases. This ef- 
fect is accounted for by introducing a density dependent 
diffusion coefficient, in terms of a simple power-law of the 
energy density D(w) ~ w"", inspired from the power-law 
nonlinearities of the model. This is the physical argu- 
ment for employing the NDE in nonlinear systems that 
exhibit some sort of chaotic diffusion. However, a rigor- 
ous derivation is yet to be found and such a calculation 
is not aimed at in this work. 
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Below, we want to present an even wider applicability 
of the NDE to 2D systems and fully resonant oscillators 
that has not been investigated before. It should be noted 
here already that this understanding of chaotic diffusion 
as the mechanism of the spreading process does not rely 
on the presence of disorder. The disorder is only essential 
for blocking non-linear waves in such system with purely 
nonlinear coupling. However, in 2D no such waves exist, 
c.f. Fig. [TJ and one can study chaotic diffusion in regular 
lattices as well, as is shown later. But first, we will an- 
alyze the 2D NDE and deduce predictions for the scaled 
spreading for two-dimensional lattices. 

The NDE in terms of a time- and space-dependent en- 
ergy density p(r, t) reads as: 



dp 
dt 



with 



pdV== E 



(5) 



Here f is the two dimensional vector and a is a nonlin- 
earity index which later will be related to the exponents 
of the spreading. The second equation represents the 
conservation of the total energy. The NDE possesses a 
radially symmetric, self-similar solution: 
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Here R{t) denotes the radius of the excitation, hence 
the excitation area of spreading states should follow 
An^ ~ R^. Although this solution only solves ([S]) for 
a source term p{f,t = 0) = E6{r) as initial condition, it 
was found that also quite arbitrary localized initial con- 
ditions asymptotically converge towards this self-similar 
solution [30|. Therefore, it can be viewed as describing 
the prototypical spreading behavior of the NDE. Substi- 
tuting B in the expression for R^, we find the scaling 
prediction of the NDE for the excitation area: 
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(8) 



with a being an unknown constant at this point. How- 
ever, for the homogeneous case k = A it is possible to 
find an analytical prediction for a as a function of the 
lattice nonlinearity k as will be presented in section ITVl 
Moreover, for a regular quadratic local potential k = 2, 
LOi^k = 1 all oscillators are in resonance and a resonant 
pertubation analysis gives an analytic prediction for a in 
dependence of the coupling nonlinearity A, also shown in 
the next section. 

For the general, nonhomogeneous (k ^ A) disordered 
case Wi^k € [0, 1] we still remain on the phenomenological 
level and check the scaling relation ([5]) by plotting the 
numerical data in appropriate coordinates. Generally, 
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FIG. 2. (Color online) Second moment An'^{t — to) for the ho- 
mogeneous nonlinearities k = A = 4 in a regular 2D (cj,;,fc — 1, 
red line) and a disordered {uji^k G [0, 1], black circle) lattice 
for /3 — 1. to was adjusted by hand to account for the long 
transient behavior before the power-law spreading is observed. 
The values for to were about to = 10*, 5 • 10* for the disor- 
dered and regular case. The dashed lines show the expected 
behavior An^ ~ — to)'^^^. Note, that here no averaging 
over initial conditions was done, the graph shows the behav- 
ior of single trajectories. In inset (a) we plot the numerical 
spreading exponent f obtained from finite difi'erences, also to- 
gether with the expectation from the NDE. Inset (b) shows 
the behavior of the structural entropy Ssti{t — to). 



one may expect that a is not a constant, but itself is a 
function of density ^/An^, in this case one would observe 
a deviation from the perfect power law ([8]) while the scal- 
ing coordinates may remain valid (such a situation was 
observed in several one-dimensional lattices [13, [2^ ) . 

From the self-similarity of the solution it immediately 
follows that the NDE predicts the structural entropy to 
be constant in the course of spreading: 



Ssti- (t) — S — S2 ^ const 



(9) 



Deviations from this relation may indicate against valid- 
ity of NDE. 



IV. NUMERICAL RESULTS 

Here, we will check the analytical predictions from the 
NDE by numerical simulations of 2D nonlinear lattices 
of the form ([T]). For the numerical time evolution we 
used a 4-th order symplectic Runge-Kutta scheme [l^l 
with step-size Ai = 0.1 {At = 0.5 in Fig. which con- 
served the energy with an accuracy of AE < 10^"^. The 
analytic scaling prediction (|5]) involves a parameter to 
that is introduced to account for transient behavior and 
is a priori unknown. We have adjusted this parameter 
in the averaged results in a way to fully emphasize the 
transient behavior of the spreading. The upper inset in 
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FIG. 3. (Color online) Second moment An^(t — to) for the 
homogeneous case k = A = 6 in a 2D lattice with disordered 
on-site potential (uJi^k £ [0, 1], black circles) and the regular 
case (uji^k = 1, red line). The transient time was adjusted 
as to — 10''. The dashed line shows the expected behavior 
An^ ~ [t — to)'^^'^ . In inset (a) we show the numerical spread- 
ing exponent u which is the slope of the curves and should 
converge to v = 3/4. Inset (b) shows the behavior of the 
structural entropy Satiit — to). Note that here no averaging 
over disorder realizations was done, the graphs show the be- 
havior of single trajectories. 

Fig.Hlexeinplarily shows the plain, unsealed results where 
no value to is adjusted. Note, that the values of to are 
several orders of magnitude smaller than the maximum 
integration times and thus do not influence the asymp- 
totic behavior at all. This parameter is solely introduced 
to increase the visibility of the scaled plots by neglecting 
the transient behavior and emphasizing the asymptotic 
spreading. 

Homogeneous nonlinearity 

At first, we chose the homogeneous case with k ~ X = 
4, 6, where the energy can be set to unity and the only 
relevant parameter is the relative coupling strength /3. 
First we note that from ([7]) we find an energy-time re- 
lation, namely t — to ^ E~"-, that is imposed by the 
scaling property of the NDE [24|. Comparing this with 
the energy-time scaling in the homogeneous (k = A) lat- 
tice we get an exact result for the nonlinear param- 
eter of the NDE in this case: a — . Hence for k = A 

Ik 

the NDE gives an exact spreading prediction, namely: 

/^n^ ^ {t-tof''/(^''-^\ (10) 

We compare this with numerical results in Figs. I2I3I 
We used a random initial excitation of 5x5 sites, where 
qi^k — and pi^k was chosen uniformly random iid., scaled 
in such a way that the total energy gave E = 1. We 
show results for (3 = 1, additionally obtained results for 



/3 = 0.1,2 with similar outcome are ommited here. For 
these simulations we show the results of a single, ex- 
emplary trajectory - hence no averaging over random 
initial conditions is performed. Contrary to the one- 
dimensional case, here the fluctuations of the propagation 
velocity are very small, because of the effective averaging 
over the large border (^ 10'^ sites) between excited and 
non-excited regions. Calculations for different random 
initial conditions with shorter integration times showed 
no significant difference. 

In both cases we find a nice convergence of the spread- 
ing towards the predicted subdiffusive law An^ ^ with 
ly = 2k/ (3k — 2), both for disordered and regular lattices. 
Insets show the convergence of the spreading exponent 
obtained from finite differences of the numerical results 
on An^. For k = A = 4, we have not quite reached the 
asymptotic regime in our simulation. The results, how- 
ever, seem to indicate a convergence towards the predic- 
tion v = 4/5. The integration for the regular lattices was 
chosen longer because of the longer transient observed 
there. Another inset shows the structural entropy which 
converges to a constant value, in consistency with the 
NDE. 

Overall, we believe that the good agreement of the nu- 
merical results with the predictions of the NDE in the 
homogeneous case k = A is a convincing evidence that 
the NDE is the correct framework to describe spreading 
in nonlinear Hamiltonian lattices also in two dimensions. 
It is quite remarkable that in this case the presence of dis- 
order does not influence the asymptotic behavior at all, 
as seen from the perfect overlap of the curves in Figs. [2] 
and[3l This shows that disorder, despite some believe, is 
not essential for chaotic diffusion. However, such a result 
can only be found in two-dimensional systems, because in 
ID disorder is required to block linear or nonlinear waves 
that would destroy the observation of diffusion. We again 
note that in this case the NDE gives an exact analytic 
prediction on the spreading exponent v which is verified 
numerically as the asymptotic behavior. 



Nonhomogeneous nonlinearity, regular lattice 

We now turn to the more general situation of k 7^ A fo- 
cusing on the case of the linear local oscillators k = 2 and 
nonlinear coupling A = 4, 6. We start with the regular 
lattice were uji^k — 1, hence all oscillators are in reso- 
nance. In this case, again an analytic prediction for the 
NDE parameter a can be obtained based on a time-scale 
analysis of the resonant dynamics. The main idea is that 
the spreading can be considered as a sequence of excita- 
tion events: given a number of already excited oscillators, 
a new oscillator on the border will be subject to a reso- 
nant forcing induced by its already excited neighbor, and 
its energy will grow due to this forcing. 

Let us consider an initially non-excited site, a neighbor 
of which demonstrates oscillations with amplitude e and 
frequency = 1 + ae^~^ (this shift of frequency follows 
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FIG. 4. (Color online) Spreading results for the inhomoge- 
neous case k = 2, A = 4 (upper graph) and A = 6 (lower 
graph) in a regular {uji^^ = 1) lattice for different energies. 
The results are averaged over random initial conditions. We 
plot the scaling prediction {An^)/E vs t/E. The upper inset 
for A = 4 shows the plain data, that is the second moment of 
M — 10 trajectories for each energy (decreasing iJ = 2 . . . 0.1) 
for lines from top to bottom in this plot ). The other insets 
show the behavior of the structural entropy Ssti{t). 



from the nonlinear coupling term). In the Hamiltonian 
for the initially non-excited site these oscillations appear 
as a driving force: 



Hi 



p2 _|_ q2 



\q — es'milt[ 
A 



(11) 



To describe the resonant excitation, we transform to the 
action-angle variables 6, I with p = —\/2I cos{9 — fit), 
q — V2Ism{9 — Qt) to obtain: 



Hi 



,A-2 7 



2/sin(6l - nt) ~esmnt\ 
A 



(12) 



Averaging over the phase of fast oscillations ilt we obtain 
a resonance averaged Hamiltonian: 



(Hi) = -ae^-^I + e^F(\/2/e~2 cos 9, \/2/e-2 sin 0) , 

(13) 



where F(a;,y) = lo^ dS, jxcos^— ysinf — sin^j'*'. The 
canonical equations of motion for /, 9 can be reduced to a 
fully dimensionless form by rescaling / e^I, t — > e^^^t, 
dropping out any e-dependence. Thus, these equations 
describe growth of energy in the driven nonlinear oscilla- 
tor to the level ~ during the time T ^ e^^^. 

Let us now consider how this characteristic time scales 
with the total energy E in the lattice. If we assume that 
the form of the distribution of the energy over the excited 
sites remains the same for all energies, then ^ E and 
we obtain that the characteristic time for the excitation 
of an initially non-excited site scales with E as: 



t - E—. 

On the other hand, for the NDE ([5]) the relation t - 
holds. Thus, we obtain: 



A- 2 



hence At 



t 



2/A 



(14) 

E-" 



(15) 



We stress here again that this analysis is based on res- 
onant excitation and definitely not possible neither for 
strongly nonlinear local potentials k > 2 nor for disor- 
dered lattices. 

We check this prediction numerically for A = 4, 6. 
We started from M = 10 (M = 24 for A = 6) differ- 
ent random intial excitations of 5 x 5 sites for each en- 
ergy {qi^k = as above), simulated the trajectories up 
a time icnd = 10^, and calculated the averages of An^ 
and S'str. Fig. |4] shows the results for the regular lat- 
tice (iOi^k = 1) for K = 2 , A = 4 (upper graph) and 
A = 6 (lower graph) for different energies. The inset in 
the upper graph shows the behavior of the second mo- 
ment for all individual trajectories, i.e. M = 10 lines for 
each energy. The asymptotically very small fluctuations 
for different initial conditions indicate that it is almost 
unneccessary to average at all for the regular case. The 
NDE predicts the energy scaling of the spreading to be 
An'^/E ~ {{t - to)/EY/^"-+^\ hence we adjust to to ac- 
count for the transient and plot Ar? jE vs (i — tQ)/E. 
The values of were around to ~ 10"^ for large energies 
and up to to = 5 • 10^ for small energies, adjusted visually 
to emphasize the power-law behavior. Indeed, this scal- 
ing gives an almost perfect collapse of data for different 
energies in both cases A = 4, 6. We also find perfect cor- 
respondance of the resonance prediction a = (A — 2)/2, 
hence An^ ^ t'' with = 2/X from above 

The structural entropy, however, shows a behavior in 
contradiction to the self-similar solution as seen in the 
insets of Fig. H) We see a clear logarithmic increase of 
S'str for all energies during the spreading. This is a clear 
indication that the spreading state is not truly self-similar 
in time. We hypothesize that this is due to the excitation 
of some long-living breather states in this system that 
trap some of the energy and decay much slower than 
the spreading process, if at all. As our measure of the 
excitation area An^ is mainly governed by the boundary 
it is insensitive to such non-decaying breathers, which 
might explain why we still find the predicted scaling. 
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FIG. 5. (Color online) Spreading results for the inhomoge- 
neous case /t = 2, A = 4 (upper graph) and A = 6 (lower 
graph) in a disordered {oJi^k £ [0, 1] iid.) lattice for different 
energies. The results are averaged over random potential re- 
alizations and the error bars indicate statistical errors. We 
plot the scaling prediction {An'^)/E vs t/E. The insets show 
the behavior of the structural entropy (5'sti )(t). 



Nonhomogeneous nonlinearity, disordered lattice 

Finally, we investigate the case of a disordered poten- 
tial oJi^k € [0, 1] also with k — 2, X = A and A = 6. In 
this case the oscillators are (typically) out of resonance 
and the predictions above do not hold. Generally, one 
expects a much slower spreading than for resonant oscil- 
lators. A conclusive analytical estimation for this case 
is still lacking at this point. Fig. [S] shows the results for 
the disordered case for k = 2 , A = 4 (upper graph) 
and A = 6 (lower graph) for different energies. Here, 
we average over 24 realizations of disorder and the total 
integration time was again tend — lO*". The disorder av- 
eraging here is essential, because of the very slow spread- 
ing which gives only about 10^ excited sites at the end 
of simulation. This is surely not enough for a reasonable 
self-averaging and hence the clarance of the results are 
improved by taken averages over disorder. The errorbars 
in Fig. U represent the statistical error of the results and 



are plotted only for the final point a.t t = 10^ for each 
energy value. The statistical error clearly increases for 
larger scaled time t/E a,s for smaller energies the states 
occupy less lattice sites and hence one finds larger fluc- 
tuations. The same argument explains why for A = 6 the 
statistical error is larger than for A = 4. Again, we follow 
the prediction on the energy scaling of the NDE and plot 
An^/E vs (t — to)/E, where again to was adjusted visu- 
ally to values around to ^ 5 ■ 10'^. As above, we find an 
almost perfect collapse of data for different energies for 
the case A = 4, for A = 6 the collapse is not as clear but 
still convincing. We observe a spreading exponent v that 
is nearly constant v^^d ~ 1/4 for A = 4 and vg^r ~ 1/6 for 
A = 6 over the studied interval, with a slight decreasing 
of ly for long times/small energy densities. 

Hence, the spreading process in the disordered case 
also obeys the energy scaling prediction of the NDE, 
moreover the saturation of the structual entropy is a sign 
of self-similarity of the spreading states. The reason why 
in regular lattices breather states emerge that destroy 
self-similarity but not in disordered lattice is still to be 
understood. One possible explanation might be that for 
a non-random potential one can apply an averaging over 
the fast oscillations of the oscillators which leads to a 
nonlinear Schrodinger lattice equation for the complex 
variable V'i.fe '■= <li,k + '-Pi.fc [ll]- This introduces a new 
conserved quantity, A/" = iV'Pj which supports the for- 
mation of breathers [2^. For a random potential, this 
averaging is not applicable, or at least leads to less exact 
results, as each oscillator has a random frequency and one 
can not apply a global averaging over these oscillations, 
hence breathers are less likely in this case. 



V. CONCLUSION 

In conclusion, we have extended the framework of 
the phcnomcnological description of energy spreading in 
strongly nonlinear lattices based on the self-similar solu- 
tion of the Nonlinear Diffusion Equation to two dimen- 
sions. This allowed for studying chaotic diffusion in the 
absence of disorder because in 2D lattice no non-linear 
waves seem to be present also for uji^k = 1- The strik- 
ing result is that one still observes chaotic diffusion very 
well described by the nonlinear diffusion equation, thus 
no disorder is required for this phenomenon. 

For the homogeneous case of equal nonlinear pow- 
ers K = A, it is possible to deduce an exact spread- 
ing law An^ ^ t^ where the exponent can be calcu- 
lated exactly from comparing the energy-time relations 
as = 2k /(3k — 2), similar to previous results in one di- 
mension [22] . Numerical simulations of two-dimensional 
nonlinear Hamiltonian lattices with k = A = 4, 6 showed 
a good convergence of the spreading towards this analytic 
expectation. 

For the non-homogeneous case, we checked the energy 
scaling prediction from the NDE for the two cases k = 2, 
A = 4, 6 for disordered and regular lattices. We found 
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that the spreading in such non-hnear lattices does show 
the predicted scahng of the number of excited oscilla- 
tors with energy and time. For the regular case there 
appears to be an additional mechanism that leads to 
deviations from the expected self-similarity of the field 
profile, as indicated by a growth of the structural en- 
tropy. For disordered potentials, however, the numerical 
results for this entropy nicely follow the scaling predic- 
tion and we conclude that in this case the NDE again 
gives the correct description of the spreading behavior. 
We have found that the effective index in the NDE for 
regular lattices (04^^ ~ 1,06, r ~ 2) is different from that 
in the disorderd case {a^^d ~ 3,a6,d ~ 5). Our explana- 
tion is that in the regular resonant mechanism 
is mostly responsible for the spreading, as neighboring 
oscillators have close frequencies. Then from the scaling 
properties of the reduced Hamiltonian which describes 
the excitation of a new site, we derived a general relation 
between the power of NDE and the nonlinearity index in 
the lattice (ITSl) which yields correct values of 04,^,06, r. 
We note that this result is not based on the assumptions 
of "strong/weak chaos" (cf. [H [ll|), but on the exact 
rescaling of the resonant Hamiltonian. For other cases 
this resonance mechanism does not work, and a theoret- 
ical derivation of the relation between a and A remains a 
challenge for future studies. 

We note that for disordered lattice in one dimension, 
a deviation from the power law behavior was found with 
an energy dependent exponent a(w), i.e. 22|. The 



results in 2D seem not to verify this claim, although a 
careful look at the results for k = 2, A = 4 in Fig. [5] 
do indicate a bending down of the curves. There are 
two reasons why such a behavior is not observed here. 
Firstly, the excitation times used in [2^ are a much better 
quantity to identify such a density dependence due to the 
different averaging, as also explained there. Moreover, 
in 2D each oscillator has not only one, but two or even 
three excited neighbors which might surpress the density 
dependence ^3^. However, the exact role of the number 
of coupled neighbors is to be examined further in later 
works. 
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